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We investigate the role of fluctuations in single molecule measurements of torque-link (t — Ik) 
curves. For semiflexible polymers of finite persistence length ( i.e. polymers with contour length 
L comparable to the persistence length Lp), the torque versus link curve in the constant torque 
(isotorque) ensemble is distinct from the one in the constant link (isolink) ensemble. Thus, one 
encounters the conceptually interesting issue of a "free energy of transition" in switching ensembles 
while making torque-link measurements. We predict the dependence on the semiflexibility parameter 
(3 — L/Lp of this extra contribution to the free energy which shows up as an area in the torque- 
link plane. This can be tested against future torque-link experiments with single biopolymers. We 
bring out the inequivalence of torque-link curves for a stiff polymer and present explicit analytical 
expressions for the distinct torque-link relations in the two ensembles and the free energy difference 
in switching ensembles in this context. The predictions of our work can be tested against single 
molecule experiments on torsionally constrained biopolymers. 
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Statistical mechanics of semiflexible polymers is of 
great current interest. Research in this area has been 
motivated by experiments [l|, Q on biopolymers in which 
single molecules are stretched and twisted to measure 
elastic properties. These experiments are designed to un- 
derstand the role of semiflexible polymer elasticity [J 0,01 
in, for instance, the packaging of these polymers in a cell 
nucleus. The process of DNA transcription can generate 
supercoihng. It is also regulated by supercoilingQ. In 
a typical experiment Q probing the twist elasticity of a 
DNA molecule, the ends of a single molecule of double 
stranded DNA are attached to a glass plate and a mag- 
netic bead. Magnetic fields are used to rotate the bead 
and magnetic field gradients to apply forces on the bead. 
By such techniques the molecule is stretched and twisted 
and the extension of the molecule is monitored by the lo- 
cation of the bead. One can thus measure the extension 
of the molecule as a result of the applied link and force 
and also make a measurement of the torque versus 
applied link |7j| . 

In the context of force-extension measurements, an iso- 
metric setup is described by the Helmholtz free energy, 
whereas an isotensional setup is described by the Gibbs 
free energy In statistical mechanics, these two ensem- 
bles are distinct [^, [l^, 11 1. In the thermodynamic limit 
these two descriptions agree, but semiflexible polymers 
(those with contour length L comparable to the persis- 
tence length Lp, i.e. f3 — L/Lp ~ 1) are not at the 
thermodynamic limit. 

In the present paper we explore this issue in the context 
of torque-link measurements. The two distinct statisti- 
cal mechanical ensembles here are the constant torque 
ensemble (isotorque ensemble) and the constant link en- 
semble (isolink ensemble). Here we focus on the role of 
fluctuations in single molecule torque- link experiments. 
In order to correctly interpret such experiments one needs 



to understand the effect of fluctuations on the measured 
quantities. For instance, it turns out, that an experi- 
ment in which the link applied to a polymer molecule 
is fixed (isolink) and the torque fluctuates yields a dif- 
ferent result from one in which the torque is held fixed 
(isotorque) and the link fluctuates [l3|. This difference 
can be traced to large fluctuations about the mean value 
of the torque or the link, depending on the experimental 
setup. Experimentally, both isolink and isotorque ensem- 
bles are realizable. Here is a schematic description of an 
experimental setup for torque-link measurements of sin- 
gle biopolymer molecules. A polymer molecule attached 
to a glass plate on one end is suspended in a suitable 
medium with a magnetized bead attached to the other 
end. The magnetized bead is kept in a magnetic trap. 
One can realize an isolink setup by using a "stiff mag- 
netic trap" and an isotorque setup by using a "soft mag- 
netic trap" which allows the applied link to fluctuate but 
applies a fixed torque to the molecule 13, 3] ■ The fiuc- 
tuations in torque-link measurements vanish only in the 
thermodynamic limit of very long polymers. In the next 
section we describe the setup in more detail. 

The paper is organized as follows. In Sec. // we dis- 
cuss the isolink and the isotorque ensembles. In Sec. /// 
we illustrate the phenomenon of inequivalence of ensem- 
bles in torque-link measurements by explicitly presenting 
analytical expressions in the context of a stiff polymer. 
In Sec. IV we draw attention to the notion of the free 
energy of transition in going from one ensemble to an- 
other and its dependence on semiflexibility which can be 
tested against future single molecule experiments on tor- 
sionally constrained polymers and simulations. Finally, 
we conclude the paper in Sec. V. 

Consider an experiment, as described in the Introduc- 
tion, in which one end of a biopolymer molecule is at- 
tached to a glass plate and the other end is attached to 
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a magnetized bead (of magnetic moment jl) kept in a 
magnetic field B which is used to rotate the bead. We 
suppose both B and fl are parallel to the glass plate. The 
energy of the bead is given by 



E = -fi.B = -fiB cos{e - 9o) 



(1) 



6*0 describes the direction of the magnetic field and 6 the 
direction of the magnetic moment fl. The variables Ik 
and Iko, which keep track of the number n of turns of 
the bead are related to 9 and 9o as follows: Ik ^ 9 + 2mT 
and Iko = + Stitt. Consider P{lk)dlk, the number 
of configurations (counted with Boltzmann weight) for a 
polymer of contour length L and bend persistence length 
Lp, characterized by a semiflexibility parameter (3 — 
in a link interval dlk of Ik. The free energy defined by 
^{Ik) = —^lnP{lk) is the free energy pertaining to a 



fixed link. The partition function [111. Il5l | for the com- 



bined system consisting of the polymer molecule and the 
magnetized bead in the magnetic field is given by: 

/ + 00 
-oo 

Constant Link Ensemble: The Limit of a Stiff Trap 
In the limit of a stiff trap {^iB oo), jl follows B 



(Ik-lkg) 



In 



closely and Ik k, Ikg. Thus cos(/fc— /fco) ~ 1 — 
this limit the partition function for the combined system 
(molccule+trap) reduces to: 

Z{lko, (3) = e^^^ d?fce'^^("=)e-'='^^^^^^ (3) 

J —oo 

Clearly, in this limit the Gaussian factor pertaining to 
the magnetic trap approaches a delta function and we 
have: 



Z(;fc, /3) « e-^*("=) 



(4) 



Here we have switched notation to write Ik in place of 
IkQ. Thus a stiff trap realizes the Constant Link (isolink) 
ensemble by constraining fluctuations 'mlk . In order to 
change the applied link from Ik to Ik + dlk one applies a 
torque < t >= Thus one gets a torque-hnk (< t > 
, Ik) curve by plotting < t > versus Ik. 

Constant Torque Ensemble: The Limit of a Soft Trap 
In the opposite limit of a soft trap /jB is small but 
large enough that the polymer does not get untwisted. In 
such a situation the link fluctuates. One can adjust Iko 
such that Ik — Iko « 7r/2. The magnitude of the torque 
\i\ = \jl X B\ = t = jiB is held fixed for a particular 
measurement. In this limit the torque t is the control 
parameter which can be changed from one reading to the 
next by changing the magnitude B of the magnetic field. 
A feedback loop is used to ensure that < Ik — Iko > is 
maintained at tt/2. Clearly, the potential energy for the 



trap, on expanding around Ik — Ikg = Tr/2 and retaining 
terms to linear order in Ik — Iko takes the form: 



E : 



-fl.B = iiB{lk ~ Iko) = t{lk - Iko) 



(5) 



Thus, in this limit Eq. ([2]) gives the following expres- 
sion for the partition function Z(i,/3) — Z(Zfco, /?)e'^*"^'' 
(where Iko is determined by the condition < Ik — Iko >= 
■|) for the combined system consisting of the polymer 
molecule and the trap 
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zit,p) 
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dlke-''^^"'^e'''"' 



(6) 



Thus in a soft trap link Ik fluctuates but torque fluctu- 
ations are constrained iJ^i]. One thus realizes the Con- 
stant Torque (isotorque) ensemble. ^(i)[l3 is the gener- 
ating function for the Ik distribution. Given the constant 
torque free energy r(i) = —^lnZ{t) one gets the mean 

link < Ik >= — ^ and the {t,< Ik >) torque-link rela- 
tion. 

Notice that Z{t) is the Laplace transform of Z{lk). 
In the thermodynamic limit of long polymers (/3 — > oo) 
the Laplace transform integral Eq.([6]) is dominated by 
the saddle point value and therefore ^{Ik) and r{t) are 
related by a Legendre transform: 



$(/fc) = T{t)+tlk. 



(7) 



For finite (3 i.e. for a polymer of finite extent, the saddle 
point approximation no longer holds true and fluctua- 
tions about the saddle point value of the free energy be- 
come important. Thus one finds that for a finite /3, $(/fc) 
and T{t) are not Legendre transforms of each other but 
are related via a Laplace transform [Eq. ©J . 

We illustrate the issue of inequivalence of ensembles 
in the context of torque-link measurements explicitly by 
considering an analytically tractable and instructive spe- 
cial case, the torque-link relation for a stiff polymer. 

Our starting point is the WLC Hamiltonian with bend 
and twist degrees of freedom in the presence of a stretch- 
ing force / and a torque t [l2| : 

H ^pI/2 + (p0 -A^)^/2 sin^ 9~ fcos9~ at"^ /2 (8) 

Here pe and are momenta conjugate to the Euler an- 
gles 9 and cj). The momentum conjugate to the Euler an- 
gle '0 is Pi/j = it, a constant of motion, which contributes 
a term — Q;t^/2 to the Hamiltonian where 

a is the ratio of the bend persistence length Lp to 
the twist persistence length Lt- The 'vector potential' 
= i<(l-cos6')[i, [13. For a stiff polymer with one 
end clamped along the z direction, we can approximate 
the sphere of directions by a tangent plane at the north 
pole of the sphere as the angular coordinate 9 always re- 
mains small. In this limit (the PWLC model [ij) where 
the tangent vector never wanders too far away from 
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the north pole of the sphere of directions, the polymer 
Hamiltoniand, 0, [H [III reduces to : 



H 



PWLC 



- A^y 

26*2 



= Hp-f- 



/(I 



-) 



where Hp is the Haniiltonian of interest in the paraxial 
limit after we take out a constant piece. In this limit 
A^, = We introduce Cartesian coordinates = 

6cos(f> and ^2 = dsm(f) on the tangent plane at the 
north pole. 

The PWLC model has been applied earlier in the con- 
text of flexible polymers at high tension Q. Here we ap- 
ply it in the context of stiff polymers at / = 0. In the 
stiff limit, the tangent vector to the polymer points es- 
sentially along a fixed direction (the north pole) even at 
f = 0. As in Ref. [4j] we restrict to the case of a = oflH. 

In Ref. [20I we have derived the partition function 
Z{f) of a stiff polymer with both end tangent vectors 
pointing along a fixed direction z at a force /. In the 
presence of a force / and a torque t this expression goes 
over to: 



P{t) = v//-tV4exp(/3/)/(27rsinh [f3^f^W/l]). 

(9) 

Notice that the "effective force" \/ {f — replaces / 

when there is a competition between a force / and a 
torque t. This is simply due to the fact that the Hamil- 
tonian Hp in the presence of a force / and a torque t 
pertains to that of a two dimensional harmonic oscillator 
with a frequency uj = yj{f — rather than lo = ^fj 

which is the expression for the frequency of a two dimen- 
sional harmonic oscillator in the presence of a force / and 
710 torque. 

At zero force, the expression simplifies and reduces to: 



Z{t) = i/(47rsin [/3i/2] 



(10) 



^From the partition function Z(t), we get the following 
analytical expression for the torque-link (t, < Ik >) rela- 
tion in the isotorque ensemble: 



(11) 



Given Z(t) one can get the partition function Z{lk) 
in the conjugate domain of an isolink ensemble (See Sec. 

//) 



Z{lk) = (7r/2/3)/(cosh [^Ik])' 



(12) 



This leads to the following torque-link (< t >,lk) rela- 
tion: 
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FIG. 1: Torque-Link Curve in the isolink (upper curve) and 
isotorque (lower curve) ensembles for (3 = 1. 



Thus, the torque-link (t, < Ik >) relation obtained in 
the isotorque ensemble (Eq.[TT]) is distinct from the one 
(< t >, Ik) obtained in the isolink ensemble (Eq.[T3]) [See 
Fig. 1]. 

In the next section we derive explicit expressions for 
the extra contribution to the free energy due to changes 
in ensemble. 

In an isotorque setup, torque is the control parame- 
ter and one measures the mean link to plot the torque- 
link {t, < Ik >) curve. In an isolink setup the roles of 
link and torque are interchanged. Consider going from a 
small torque {lki,ti) to a large torque (^^2,^2) configu- 
ration via an isotorque setup and returning from a large 
link (lk2, 12) to a small link {lki,ti) configuration via an 
isolink setup. Since the torque-link relation depends on 
the chosen ensemble, in general there will be two dis- 
tinct curves joining the points {lki,ti) and [Ik^^t^) in 
the torque-link plane, describing the two processes. Such 
a transformation could lead to a net area being enclosed 
in the torque-link plane. This appears paradoxical since 
it seems to suggest that one can extract work from the 
system via a cyclic process. The resolution of this para- 
dox is as follows. In completing the cycle and returning 
to the initial state one is in fact changing ensembles twice 
at the two end points. Such ensemble changes in a cyclic 
transformation involve finite changes in free energy which 
need to be taken into consideration. In particular, we no- 
tice for the special case of a stiff polymer the torque-link 
relations (Eq. [llj and Eq. |13| ) lead to a free energy dif- 
ference of : As4i = 



2 cosh(7r^fc2) | , N , 1 

TTTTT ~(*2ifc2-ii/fci)-t-- 

p cosh(7r«fci) p 



'In 



sin(g^) 
sin(^) 



< t >= (27r//3) tanh [-Klk]. 



(13) 



for a transformation between the states {lki,ti) and 
(^^2,^2)- In the stiff regime (i.e. at small /3) the de- 
pendence of the free energy on /? will be dominated by 
the first term. In other words, the free energy difference 
^stiff ^ ^- This is a prediction of our analysis which 
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can be tested against experiments with stiff biopolymers 
like actin filaments. 

An analysis similar to the one in Ref. 0] applied to the 
context of torque-link measurement shows that one gets 
a contribution A = j^ln^ {Ik^,} to the free energy com- 
ing from fluctuations around the long polymer (/? — > oo) 
limit by expanding A{lk) = — tlk around the sad- 

dle point value Ik = Ik^, pertaining to the long polymer 
limit. Here $(^fc) corresponds to the isolink free energy. 
This extra contribution A{lk) vanishes in the limit of 
(3 — > oo. For finite /3, this nonzero contribution to the 
free energy accounts for the transition between the con- 
stant link ensemble and the constant torque ensemble. 
Work is done on the bead by the trap in making the 
trap stiffer, while in going from a stiff trap to a soft trap 
work is extracted from the bead by the trap. The net 
work done is the difference between the work done at 
the two ends of the torque-link curves in switching en- 
sembles. This net work exactly cancels out the nonzero 
area enclosed in the torque-link plane. The net area en- 
closed in the torque-link plane pertaining to the "Free 
Energy of Transition" scales as 1//3 and therefore grows 
with the rigidity of the polymer. These predictions of 
our study can be tested against future simulations and 
single molecule experiments. 

To summarize, in this paper we have studied inequiv- 
alence of ensembles for torque-link measurements. We 
have calculated the free energy difference between torque- 
link measurements in the isotorque and isolink ensembles 
for a stiff polymer. In addition, we have determined the 
contribution to the "free energy of transition" in going 
between the isolink and the isotorque ensembles by ex- 
panding the free energy difference A{lk) around the long 
polymer limit. We predict the dependence on the semi- 
flexibility parameter P of this extra contribution to the 
free energy which shows up as an area in the torque- 
link plane. For the special case of a stiff polymer we 
find explicit analytical expressions for the torque-link 
(t, < Ik >) relation obtained in the isotorque ensem- 
ble (Eq.[ll]) and show that it is distinct from the one 
(< t >,lk) obtained in the isolink ensemble (Eq.[T3]). We 
also show that in this stiff regime the free energy differ- 
ence Astiff ~ ^- AH the predictions mentioned here can 
be qualitatively and quantitatively tested against future 
single molecule experiments on torsionally constrained 
biopolymers. 

The theoretical predictions presented in this study are 
expected to generate interest in torsionally constrained 
single molecule experiments which will eventually lead 



to a deeper understanding of the role of twist elasticity 
in biological processes involving gene regulation 22 1. 
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